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Abstract. 

In  this  paper,  the  propagation  of  interfacial  waves  near  the  critical  depth  level  in  a 
two-layered  fluid  system  is  investigated.  We  first  derive  an  evolution  equation  for  weakly 
nonlinear  and  dispersive  interfacial  waves  propagating  predominantly  in  the  longitudinal 
direction  of  a  slowly  rotating  channel  with  gradually  varying  topography  and  sidewalls. 
The  new  evolution  equation  includes  both  quadratic  and  cubic  nonlinearities.  For 
interfacial  waves  propagating  in  certain  type  of  non-rotating  channels  with  varying 
topography,  we  find  two  families  of  periodic  solutions,  expressed  in  terms  of  the  snoidal 
function,  to  the  variable-coefficient  equation.  As  the  limiting  cases  of  these  periodic-wave 
solutions,  a  family  of  solitary- wave  solutions  and  an  isolated  shock-like  wave  solution 
are  also  obtained.  In  a  uniform  rotating  channel,  our  small-time  asymptotic  analysis 
and  numerical  study  show  that  depending  on  the  relative  importance  of  the  cubic 
nonlinearity  to  quadratic  nonlinearity,  the  wavefront  of  a  Kelvin  solitary  wave  will 
curve  either  forwards  or  backwards,  trailed  by  a  small  train  of  Poincare  waves.  When 
these  two  nonlinearities  almost  balance  each  other,  the  wavefront  becomes  almost  / 
straight-crested  across  the  channel,  and  the  trailing  Poincare  waves  diminish. 
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1.  Introduction 

Nonlinear  and  dispersive  internal  wave  trains  have  been  observed  and  measured 
in  marine  straits  (e.g.  Gargett  1976;  Alpers  &  Salusti  1983;  La  Violette  &:  Arnone 
1988)  and  in  long,  narrow  thermally  stratified  lakes  (e.g.  Thorpe,  Hall  &  Crofts  1972; 
Hunkins  &  Fliegel  1973;  Farmer  1978).  Many  attempts  have  been  made  to  extend 
the  applicability  of  the  original  Kadomtsev-Petviashvili  (KP)  equation  (Kadomtsev  & 
Petviashvili  1970)  to  accommodate  important  factors  affecting  these  observed  internal 
wave  propagations,  such  as  rotation,  variations  of  topography  and  sidewalls,  and  a 
background  current  field.  Most  recently,  Chen  &  Liu  (1995)  presented  the  derivation 
of  the  unified  KP  (uKP)  equation  for  weakly  nonlinear  and  dispersive  interfacial  waves 
of  two-layer  fluids  propagating  predominantly  in  the  longitudinal  direction  of  a  slowly 
rotating  channel  with  gradually  varying  topography  and  sidewalls,  and  a  weak  steady 
background  current  field.  Their  result  was  a  generalization  of  all  previous  work  within 
the  context  of  surface  and  interfacial  waves. 

The  uKP  equation  was  derived  under  the  assumption  that  the  effects  of  nonlinearity,' 
dispersion,  rotation,  transverse  modulation  and  variations  of  topography  and  sidewalls 
are  weak  but  equally  important.  In  the  equation,  the  nonlinear  term  is  quadratic  and 
its  coefficient  is  proportional  to  D_2  =  p~ l{h~y  —  l{h^y,  where  p^  and  are 
the  dimensionless  densities  and  leading-order  depths  in  the  upper  and  lower  layers, 
respectively.  This  implicitly  limits  the  range  of  the  variation  of  topography  because  D-2 
must  remain  same  sign  and  be  an  order  one  quantity  to  ensure  that  the  nonlinearity 
indeed  is  as  important  as  the  dispersion  and  other  effects.  If  the  bottom  changes 
to  certain  extent  that  £>_2  changes  sign,  the  uKP  equation  is  no  longer  applicable 
because  near  the  critical  depth  level  defined  as  =  (p'*'//!>~)^/^(i.e.  Z)_2  =  0),  the 

coefficient  of  the  nonlinear  term  is  so  small  that  the  balance  between  the  nonlinearity 
and  dispersion  (and  other  effects)  becomes  impossible  under  the  Boussinesq  assumption, 
i.e.  O(ao/fio)  =  f^(^o/^o)  ^  where  uq,  and  /o  are  the  typical  wave  amplitude,  depth 


and  wavelength,  respectively. 

In  the  study  of  the  change  of  polarity  of  two-dimensional  solitary  waves  as  they 
pass  through  the  critical  depth  level,  Helfrich,  Melville  &  Miles  (1984)  showed  that  near 
the  critical  depth  level,  nonlinearity  and  dispersion  can  be  of  the  same  order  magnitude 
in  the  parametric  regime  where  0{aofho)  =  0{hollo)  ^  1.  In  this  situation,  the  cubic 
nonlinearity  becomes  comparable  to  or  dominates  the  quadratic  nonlinearity  and  must 
be  taken  into  account. 

The  inclusion  of  the  cubic  nonlinearity  can  have  a  significant  influence  on  wave 
motion.  For  transcritical  two-layer  flow  over  a  two-dimensional  topography  or  past  a 
constriction  in  a  stationary  channel,  Melville  and  his  associates  found  that  the  inclusion 
of  the  cubic  nonlinearity  changes  the  character  of  the  upstream  disturbance,  giving  an 
upstream  advancing  monotonic  non-dissipative  bore,  rather  than  a  trains  of  Boussinesq 
solitary  waves  (Melville  &;  Helfrich  1987;  Tomasson  &;  Melville  1991). 

Although  recently  there  have  been  several  studies  on  Kelvin  solitary  wave 
propagation  in  a  rotating  channel  (Katsis  &  Akylas  1987;  Melville,  Tomasson  &  / 

Renouard  1989;  Grimshaw  &  Tang  1990;  Chen  &  Liu  19966),  Kelvin  solitary  wave 
propagation  near  the  critical  depth  level  has  not  been  studied  yet.  We  expect  that 
the  inclusion  of  the  cubic  nonlinearity  may  give  rise  to  some  new  wave  propagation 
phenomena,  especially  when  the  Kelvin  solitary  wave  propagates  over  a  vaxying 
topography  such  that  the  coefficient  of  the  quadratic  nonlinear  term  changes  its  sign. 

In  this  paper,  we  consider  the  formulation  and  analytical  solutions  of  the  evolution 
equation  that  governing  three-dimensional  interfacial  wave  propagation  near  the  critical 
depth  level  in  the  most  possible  general  setting.  Then,  we  apply  the  equation  to  study 
Kelvin  solitary  wave  propagation  near  the  critical  depth  level  in  a  rotating  uniform 
channel  to  reveal  new  wave  phenomena  caused  by  the  cubic  nonlinearity. 

In  the  next  section,  assuming  that  the  effects  of  nonlinearity,  dispersion,  rotation, 
transverse  modulation  and  variations  of  topography  and  sidewalls  are  weak  but 
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equally  important,  we  derive  an  interfacial  displacement  evolution  equation  for  waves 
propagating  near  the  critical  depth  level.  In  section  3,  we  study  the  integrability  of 
the  evolution  equation  and  seek  analytical  solutions  to  describe  wave  propagations  in 
certain  type  of  curved  channels  with  varying  topography.  In  section  4,  we  carry  out 
the  small-time  asymptotic  analysis  and  numerical  study  of  the  evolution  equation  for 
Kelvin  solitary  wave  propagation  in  a  rotating  uniform  channel  to  examine  the  physical 
features  brought  by  the  cubic  nonlinearity.  Concluding  remarks  are  given  in  section  5. 


2.  Derivation  of  the  evolution  equation 

We  consider  interfacial  waves  propagating  along  the  interface  of  two  fluid  layers 
confined  in  a  rotating  wide  channel  with  varying  topography  and  sidewalls.  Cartesian 
coordinates  are  employed  with  the  x-  and  y-axis  in  the  longitudinal  and  transverse 
directions  of  the  channel,  respectively,  and  the  2-axis  pointing  vertically  upwards.  The 
interface  is  denoted  by  2  =  fj{x,y,t)  with  2  =  0  being  the  undisturbed  interface.  The 
rotation  rate  about  the  vertical  axis  is  f/2.  The  densities  of  the  upper  and  lower  layers 
are  p"*"  and  p~  (p~  >  p"*"),  respectively  (hereafter  superscripts  -f  and  —  are  used  to 
identify  quantities  in  the  upper  and  lower  layers,  respectively).  The  upper  and  lower 
layers  are  bounded  by  2  =  (the  rigid-lid  assumption  is  adopted  to  approximate  the 
free  surface)  and  2  =  H~(x,y),  respectively.  The  vertical  sidewalls  are  represented  by 
y  =  yi{x)  and  y  =  yr{x)  with  yi  >  yr- 

The  fluids  in  both  layers  are  assumed  to  be  inviscid  and  incompressible.  The 
dimensionless  governing  equations  and  boundary  conditions  for  flows  in  the  upper  and 
lower  layers  are 

(2.1a) 


du^  dv^  dw^ 

+  —  +  —  =  0, 


dx  dy 


dz 


du^ 

dt 


+  e  « 


.du^ 

dx 


-l-u" 


.du^ 

dy 


-|-  Xif 


.du^' 

dz  , 


—  7U 


dp^ 


(2.1b) 
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dv^  (  ,  dv^  .  dv^  .  dv^  \  . 

+  e  1  u^— - hu  -7;—  +  t«  1  +7^^  = 


dt 


dx 


dw^ 

dt 


+  e\u 


.dw^ 

dx 


dy 


dz 


.dw^ 

dy 


+  w 


.dw^ 


dp^ 
dy  '' 

dp^ 

'~d7' 


+  dy  {  +dri  +dr]^ 

^  =Tt+Ti*^ry, 


on  2  =  e?/, 


—  p  p  +  77  =  0  on  z  =  ey^ 
w'^  =  0  on  2  =  ^ 


w  =  —u 


dH- 

dx 


dH- 

dy 


on  z^-H  (x,?/), 


on  y  -  yr{x),yi{x), 


(2.1c) 

(2.1d) 

(2.1e) 

(2.1f) 

(2.1g) 

(2.1h) 

(2.1i) 


where  superscripts  ±  are  vertically  ordered;  and  are  the  dimensionless  velocity 

components  in  the  x-,  y-  and  2-direction,  respectively,  and  p*  is  the  dimensionless 
hydrodynamic  pressure.  The  dimensionless  variables  are  related  to  their  dimensional 
counterparts  (denoted  by  tildes)  by: 


(x,y)  =  lo{3:,y),  z  =  hoz,  i  =  —t,  =  hoH^,  Pr  =  kyr, 

Co 

-  7  ~  <^0^3/  ±  ±S  O^oQ)  ± 

yi^loVh  V  =  ^or],  [U  ^v^)  =  -—{u  ,v  ),  w  = -j—w  , 

riQ  Iq 


pop  ,  Co  =  {p  -p+)p/jo,  P 


^0  2  ±  ± 
T^oPoP  P  , 

flo 


(2.2) 


where  Iq  and  ho  are  the  characteristic  wavelength  and  depth,  respectively;  oq  and  cq  are 
the  characteristic  amplitude  and  phase  speed  of  linear  long  waves,  respectively;  g  is  the 
gravitational  acceleration  and  po  is  the  characteristic  density.  The  parameters  e,  and 
7  appearing  in  (2.1)  are  defined  as 


(  =  ao/ho,  P  =  (ho/l<iY,  7  =  <o//«o. 


(2.3) 


which  measure  nonlinearity,  dispersion  and  rotation,  respectiveh^ 

We  shall  derive  the  evolution  equation  for  weakly  nonlinear  and  dispersive  interfacial 
waves  propagating  near  the  critical  depth  level  and  travelling  predominantly  in  the 
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longitudinal  direction  of  a  slowly  rotating  channel  with  gradually  varying  topography 
and  sidewalls.  More  specifically,  we  assume  that 


D-2  =  P  /{h  y  -  p'^ /{h^y  =  0(e),  e  =  ap,  7  =  fip,  with  p  <  1,  (2.4) 

where  a  and  (3  are  two  order  one  constants,  and  the  bottom  and  the  sidewalls  can  be 

expressed  as 

H~  =  h~{p'^x)  +  (2.5a) 

Vr  =  -Yr{p'^x),  yi  =  -Yl{p^x).  (2.5b) 

fX  fi 

We  use  p  as  the  basic  perturbation  parameter  in  the  following  derivation. 

To  derive  a  single  evolution  equation  for  the  interfacial  displacement,  we  introduce 
the  following  transformation; 

^  C~^{p^x)dx-t,  X  =  p^x,  Y  =  py,  Z  =  z,  (2.6a) 

where  C  is  given  by 

C{X)  =  (p+/fi+  +  p-jh-)~^'\  (2.6b/ 

and  h'^  =  H'^  =  constant  is  used  for  convenience.  Note  that  C  is  the  leading  order  of 
the  local  linear-long-wave  speed  and  ^  =  0(1)  is  the  characteristic  coordinate  moving  at 
the  speed  of  C  to  the  right. 

In  the  moving  coordinates  (^,  X,Y,Z),  governing  equations  and  boundary  conditions 
(2.1)  can  be  rewritten  as 


Y  =  py. 

Z  =  z, 

(2.6a) 

(2.6b/ 

1  du^  ^  n 


(2.7a) 


du^ 

_L 

.(\du^ 

a[x  u 

Ic  a?  + 

1  dp^ 

C  di 

ax' 

dv^ 

-t 

( 1  dv'^ 

dt 

“  afi  u 

1 

dXj 

+  aF  +  "’  az 


+  dv^  .  dv^ 

ay +“■  az 


I3^u^  =  fi- 


(2.7b) 

(2.7c) 
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2  f 


u 


'  1  dw^ 


,dw^' 

dX 


+  fj.v^ 


.dw^ 

'W 


+  w^ 


.dw^ 


dZ 


dp^ 

1z 


± 

w  =  -—  +  afi 


+  f  I  dp  2dp\  ±dp 

“  (caf+'‘ax)+'‘’’V 


on  Z  =  e?7, 


w 


-pZ’u' 


p^p'^  —  p  p  +  ?7  =  0  on  Z  =  ep^ 


dB 

dY 


«;+  =0  on 

dfi- 

,  2dB\  3  . 

dx 

± 

^dF 

V  : 

=  pu  —  on 

(2.7d) 

(2.7e) 

(2.7f) 

(2.7g) 

(2.7h) 

(2.7i) 


where  (2.5a)  and  (2.5b)  have  been  used. 

A  solution  to  the  governing  equations  and  boundary  conditions  (2.7)  is  sought  in 
the  following  series  forms: 

G  =  Go(e,  F,  Z)  +  /xGx(^,  X,  F,  Z)  +  A,  F,  Z)  +  0{p% 


A,  F,  Z)  +  A,  F,  Z)  +  0{p% 

p  =  77o(e,  A,  F)  +  pp,  (e,  A,  F)  +  p-^p^ii.  A,  F)  +  0{p% 


(2.8a) 

(2.8b)^ 
(2.8c) 

where  G  =  Substituting  (2.8)  into  (2.7)  and  expanding  the  interfacial 

boundary  conditions  (2.7e)  and  (2.7f)  at  Z  =  0  and  the  bottom  boundary  condition 
(2.7h)  at  Z  =  — h",  we  obtain  a  sequence  of  initial-boundaxy-value  problems  by 
collecting  coefficients  of  ju".  Carrying  out  the  perturbation  analysis  to  the  second  order 
(n  =  2),  we  find  out  that  the  compatibility  condition  of  the  second-order  problem 
requires  po  and  pi  to  satisfy  the  following  evolution  equation  (for  details,  see  Chen 
1995): 


-1/2  d  f^indpo]  3aG2  ^  d'^pf 

^  W) 


—  _  2p2  j-\  d^pQ  Di  d'^po 

D-2  0^2  ^  ^  -^-3  fi(2  +  g 


C^d^Po  (3^  p-BC'^d'^po  3aC^  ddpopi  ^ 

^  2  5F2  2^^°'^  2(h-)2  5^2  2  a^2 


2 


(2.9) 
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where  Di,  D_2  and  D-z  are  defined  as 


Dr.{x)  =  p-{h-Y  +  (-ir v+(h+r  (n  =  1,  -2,  -3). 


(2.10) 


Because  Z)_2  =  p~ l{h~Y‘  —  p'^ =  0(e),  the  last  term  involving  rji  in  (2.9)  can 
be  dropped  and  we  obtain  the  evolution  equation,  including  both  quadratic  and  cubic 
nonlinearity,  for  interfacial  waves  propagating  near  the  critical  depth  level: 


pm  ^  ( ^172^^0 

dx  1  di 


7^  ,  p-BC^d^iio  „ 


0^5%  7^  p  BC^d^ 
2  dY^  2(h-)2  ap 


(2.11) 


The  leading  order  of  the  impermeable  boundary  conditions  along  the  sidewalls, 
(2.7i),  provides  the  sidewall  boundary  conditions  for  ‘i]q: 


drjo  7  _  1  dF57o 

dY^  CdXdi 


onY  =  Yr{X\Yl{X). 


(2.12) 


Note  that  if  0(e)  =  0{p?)  and  Z)_2  =  0(1),  the  cubic  nonlinear  term  in  (2.11) 
becomes  higher-order  and  can  be  neglected;  the  resulting  equation  recovers  the  uKP 
equation  derived  by  Chen  &  Liu  (1995)  through  the  following  changes: 


(2.13) 


The  boundary  conditions  on  the  sidewalls  for  the  uKP  equation  are  also  recovered 
from  (2.12)  through  the  same  changes.  Besides  the  cubic  nonlinear  term,  the  physical 
meaning  of  each  of  the  remaining  terms  in- (2.11)  is  the  same  as  that  of  the  corresponding 
term  in  the  uKP  equation  (Chen  &  Liu  1995). 

In  the  absence  of  rotation,  (2.11)  can  be  further  reduced  to  the  existing  evolution 
equations  for  interfacial  wave  propagation  near  the  critical  depth  level.  For  example,  if 
the  bottom  is  flat,  (2.11)  becomes: 


^2^3  Did^r]Q 

dXdi'^  "^5^  p'^^  de  ^  2  dY^ 


(2.14) 
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which  agrees  with  the  modified  KP  equation  given  by  Tsuji  &  Oikawa  (1993), 
who  numerically  solved  the  equation  to  study  the  interaction  of  two  solitary  waves 
propagating  in  close  directions.  On  the  other  hand,  if  the  wave  field  is  independent  of  Y 
(the  sidewalls  should  be  straight  and  parallel),  (2.11)  reduces  to 


3e 


di 


D.d^Vo  p-BC^dr,o_ 

6  de  2(/i-)2 


When  .6  =  0,  (2.15)  is  equivalent  to  the  evolution  equation  derived  by  Helfrich  et  al. 
(1984).  When  h~  is  constant,  (2.15)  is  identical  to  the  evolution  equation  derived  by 
Chen  &  Liu  (1996a)  to  study  the  effect  of  a  random  topography  on  two-dimensional 
interfacial  wave  propagation. 

Assuming  that  the  solution  of  (2.11)  with  boundary  conditions  (2.12),  t/q,  is  locally 
confined,  we  find  out  that 


X  =  {Ciwf^  r"  T^od^dr,  W{X)  =  Yl{X)  -  Yr{x)  (2.16) 

JYr  J—OO 

is  the  first-order  invariant  when  the  rotation  is  absent  (i.e.  7  =  0)  and  / 

[Yl  r+oo 

J  =  C  /  rild^dY  (2.17) 

JYjl  J-~-oo 


is  the  second-order  invariant  regardless  whether  the  rotation  is  present.  The  locally 
confined  solution  also  has  to  satisfy  the  same  constraint  imposed  on  the  uKP  equation 
(Chen  &  Liu  1995): 

/+00 

Tjod^  =  F{X)  exY>{-l3YIC).  (2.18) 

-OO 

Thus,  again,  the  minimal  requirement  for  a  solution  of  (2.11)  with  (2.12)  to  be  locally 
confined  is  that  its  initial  condition  satisfies  (2.18)  (with  /?  =  7/p)  at  X  =  0. 


3.  Analytic  solutions 

Equation  (2.11)  is  a  nonlinear  evolution  equation  with  variable  coefficients.  In 
general,  no  procedure  is  available  to  obtain  analytic  solutions  of  the  corresponding 
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Cauchy  problem,  let  alone  initial-boundary-value  problems.  In  this  section,  we  focus  on 


initial- value  problems  first. 

A  powerful  tool  to  investigate  the  complete  integrability  of  a  nonlinear  evolution 
equation  is  the  Painleve  PDE  test  (Ablowitz  &  Clarkson  1991).  It  provides  a  useful 
criterion  for  whether  a  given  partial  differential  equation  is  completely  integrable.  We 
now  carry  out  the  Painleve  analysis  to  find  out  under  what  conditions  (2.11)  passes  the 
test,  as  we  did  for  the  uKP  equation  in  Chen  &  Liu  (1995). 

To  simplify  the  algebraic  manipulation,  we  introduce  the  following  transformation: 


iE  ( 

tC  UP-s 


c,  .  =  / 


which  transforms  (2.11)  into 


AL 

drd( 

where 


+  -  2^  +  P  +  +  Hr.  Y)^  +  c(r)^  -  d(r)(  =  0,  (3.2) 


a=:3CyDu  b  = 


Zp-BC^ 


zp-c  \{h-Y  3  c^h-yi  dh- 
Di{h-y[  Di  2  dX’ 


ZCD.2 


2p  \DiD-3 


(3.3a)_ 


(3.3b) 


(3.3c) 


We  seek  solutions  of  (3.2)  in  the  form 

OO 

r{T,i,Y)  =  (  +  i,{T,Y),  (3.4) 

i=0 

where  p  is  an  integer,  and  Y)  and  Uj(r,  Y){j  =  0, 1, 2, . . .)  are  analytic  functions  of  r 
and  Y  in  the  neighborhood  of  a  non-characteristic  movable  singularity  manifold  defined 
by  (^  =  0.  Upon  substitution  of  (3.4)  into  (3.2),  the  leading-order  analysis  requires  that 
p  =  —1  and  Mg  =  1.  Equating  the  like  powers  of  tp  yields  the  general  recursion  relation: 


(3.5a) 


where 


W,  = 


a^-JkA  -  duj-4  -  2(j  -  3)(j  -  4) 


dV^ 


j-1  k  i-1 

k=l  m=:0  m=l 


i-1 

■E 

k^O 


+  e(i  -  3)(i  -  4)  +  (j  -  3)(i  -  4)uj. 

+  (i  -  4) 


d^b  (di)V  ' 


/  ,  ,  d^dUj^z 

“.-3  c  +  Ojpj  1  + -^  +  2<. 


(3.5b) 


for  j  >1  (define  Uj  =  0  for  j  <  0). 

Since  J  =  4  is  a  double  root  of  the  indicia!  polynomial  determining  the  resonances 
(occurring  at  some  j  where  Uj  is  arbitrary),  (3.2)  cannot  pass  the  Painleve  test.  Note 
that  the  appearance  of  the  double  root  is  due  to  the  existence  of  the  cubic  nonlinear 
term  in  (3.2);  even  for  the  simplest  case  when  the  rotation  is  absent  and  the  bottom  is 
fiat,  the  equation  still  cannot  pass  the  Painleve  test.  Since  there  is  no  way  to  let  (3.2) 
pass  the  Painleve  test  by  imposing  restrictions  on  the  variable  coefficients  a,  5,  c,  d  and 
e  (if  there  were,  the  evolution  equation  could  be  reduced  to  a  completely  integrable 
equation  in  canonical  form  via  an  elementary  transformation),  the  search  for  analytical  -^ 
solutions  of  (3.2)  or  (2.11)  becomes  very  challenging.  In  what  follows,  we  shall  give  some 
analytical  solutions  to  (2.11)  for  certain  type  of  varying  topography  and  sidewalls  in  the 
absence  of  rotation. 

Under  the  assumption  that 


/?  =  0,  h  =  const. 


o~  BC 

■^  =  F,(X)Y  +  F„{X), 


2(h-f 

where  Fq  and  Fi  are  arbitrary  functions  of  X,  (2.11)  becomes: 


(3.6) 


dXdi 


-a^CD.z^  + 


Di  d% 


5^2  '  6C 

+  f5^  +  [F.(A')l'+FoW]^  =  0, 


(3.7) 


which  is  still  a  nonlinear  evolution  equation  with  variable  coefficients.  Assumption  (3.6) 
implies  that  there  is  no  rotation  and  the  variation  of  topography  is  weak  and  behaves 
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like  a  linear  function  in  the  transverse  direction.  Note  that  this  assumption  actually  is 
the  complete  integrability  conditions  for  the  uKP  equation  (Chen  &  Liu  1995). 

We  find  that  the  following  transformation: 


rjo  = 


1 

aC 


(3.8a) 


f  =  Vqd,x/c, 
rX  r  C 


X  =  ^/6  -  r  jT  Fr{q)dq  -  |Fo(9)  +  |  Fi(s)ds)  J  d<?| . 

Y  =  YjC  +  Fy[s)dsdq  , 

transforms  equation  (3.7)  and  boundary  conditions  (2.12)  with  7  =  0  (i.e.  (d 

dX  Vaf  ^  ^dX  ^  dX  ^  dX^ )  ^  ^  dY^ 


with 


and 


CD.2 


r  = 


^  =  (6A)-'^ 


dW 


R 


dX 

dYL 

dX 


+  C 


d-C 


r  F,{q)dq 
Jo 

I 

Jo 


dC 


dX 

K 

dX 


^  on  Y^Yn{Tl 


on  y  =  yh(f), 


respectively,  where 


Yr{T)  = 

YL{f)  =  dM 


yR/(^  +  l  l'F,{s)dsdq 

Yl/C  +  £  F,{s)dsdq 


(3.8b) 

(3.8c) 

(3.8d) 
0)  into 

(3.9) 

(3.10) 

f 

(3.11a) 

(3.11b) 


(3.12a) 

(3.12b) 


As  shown  in  (3.10),  the  parameter  r  signifies  the  relative  importance  of  the  quadratic  and 
cubic  nonlinearities.  Note  that  with  a  rescaling,  (3.8)  is  equivalent  to  the  transformation 
used  to  transform  the  uKP  equation  (under  assumption  (3.6))  into  the  KP  equation  in 
canonical  form  (Chen  h  Liu  1995). 


Equation  (3.9)  is  a  constant-coefficient  modified  KP  equation.  Although  the  KP 
equation  and  the  modified  KdV  equation  are  both  completely  integrable,  the  modified 
KP  equation  (3.9)  is  not  (note  that  (3.9)  is  a  special  case  of  (3.2)).  Neverth^ess, 
since  (3.9)  is  a  constant-coefficient  wave  equation,  it  may  admit  travelling-wave 
solutions.  We  now  seek  for  all  bounded  travelling- wave  solutions  of  (3.9).  Substituting 
Ci4>)  =  +  lY  —  ojf)  (I  and  u  are  constants)  into  (3.9)  and  integrating  the  resulting 

equation  twice,  we  obtain 


C"  =  -  3rC"  +  {u;-  3/2)C  +  A6  +  B.  (3.13) 


where  '  =  d/d</),  A  and  B  are  two  constants  of  integration.  For  bounded  solutions,  A 
must  be  zero.  Therefore,  (3.13)  can  be  recast  into 

X'  =  ‘d,  ■&'  =  X^-PXAq,  (3.14a) 

where 

C  =  X  +  ^5  P  =  +  3/^  —  a;,  g  B  +  (w  —  3/^)r  —  2r^.  (3.14b) 

> 

Carrying  out  the  analysis  in  the  phase  plane  (x,^)  (for  details,  see  appendix),  we  find  ' 
two  families  of  periodic  solutions  (corresponding  to  close  orbits  in  the  (y,  i?)-plane)  for 
(3.9),  which  are  given  as  follows. 

The  first  family  of  periodic  solutions  to  (3.9)  is  given  by 


with 


C  =  c  —  i/[m(3  —  2n)  —  n(2  —  n)][4n(m  —  n)]  ^ 

-H  i?(l/n  — 1)  l-nsn^($,m)]  \  (3.15a) 

$  =  K  sec  -f  $0  =  A'  sec  Oq  ^cos  OqX  d-  sin  OqY  —  VT ^  +  $o,  (3.15b) 

K  =  H^J{1  -  n)[8n(m  -  n)]-i,  (3.15c) 

V  =  cos  6q  |3r^  +  3  tan^  —  i?^[m^(3  —  4n  -|-  4n^) 

—  2mn(2  —  n  -|-  2n^)  -t-  n^(4  —  4n  -|-  3n^)][4n(m  —  n)]“^| , 


(3.15d) 
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where  H  is  the  wave  height,  m  and  n  (0  <  n  <  m  <  1)  are  two  parameters,  6o  is  the 
angle  of  incidence,  K  is  the  X-component  wavenumber,  V  is  the  phase  speed  and  #o  is 
a  phase  constant.  The  wavelength  of  the  periodic  wave  (3.15),  A,  is  given  by 


A  =  2K(m)  cos  6o  /  K, 


where 


K(m)  =  f 

Jo 


(3.16) 


(3.17) 


'  '  Jo  [l-msmn)>/2 

is  the  complete  elliptic  integral  of  the  first  kind. 

Assuming  that  the  mean  interface  is  located  at  ^  =  0,  i.e.  the  net  area  covered  by 


the  interface  within  a  wavelength  is  zero: 


/2K(m)  -  A  /•2K(«) 


(3.18) 


we  obtain 


(1  —  l/n)n(n,  m)/K(m)  =  r/i/ —  [m(3  —  2n)  —  n(2  —  n)]  [4n(m  —  n)]  (3.19) 


where 


p7l 

n(n,  m)  =  j 
Jo 


(3.20) 


’  Jo  [1  —  nsin^  t][l  —  m  sin^  ''  ' 

is  the  complete  elliptic  integral  of  the  third  kind. 

Upon  using  (3.16)  and  (3.19),  (3.15a)  can  be  rewritten  as 

C  =  ^(1  -  1/w)  |n(n,m)/K(m)  -  1  -  nsn^  |.  (3.21) 

Given  H  (wave  height),  A  (wavelength)  and  (angle  of  incidence),  we  first  express  n  in 
terms  of  m  through  (3.15c)  and  (3.16): 


i  [(m  +  r)  ±  ^(m  +  r)2-4r]  , 


(3.22a) 


where 


(3.22b) 


Note  that  since  n  (0  <  ?2  <  1)  is  a  real  parameter,  the  range  of  m  is  a  subset  of  (0,1)  such 
that  the  determinant  [(m  +  F)^  —  dF]  >  0.  Substituting  (3.22)  into  .(3.19),  we  obtain  a 
transcendental  equation  for  m.  It  can  be  showed  that  for  r  <  0(r  >  0),  no  solution  can 
be  found  when  n  =  n+(n  =  7r_)  is  substituted  into  (3.19),  whereas  a  unique  solution 
is  obtained  when  n  =  n_(n  =  n+)  is  substituted.  Indeed,  for  the  same  value  of  m,  if 
{m,n-)  corresponds  to  (C, f"),  then  should  correspond  to  (— C,—’')  with  a  phase 

shift  K(77r).  The  reason  is  that  (3.9)  is  invariant  under  transformation  (C,  ?*)  ^  ("Ci  ~^)- 


Figure  1.  Profiles  of  the  periodic  wave  (3.21)  for  different  values  of  Pi  with  P2  =  —0.5: 

( - ),  Pi  =  5.0  (m  =  0.7229,  n  =  0.2869);  ( - ),  Pi  =  10.0  (m  =  0.9549,  n  =  0.3856); 

( - ),  Pi  =  15.0  (m  =  0.9943, n  =  0.4507).  For  comparison,  the  sinusoidal  wave  profile 

(3.24)  is  also  plotted  in  this  figure  ( - ). 

From  (3.19)  and  (3.22),  one  can  see  that  the  normalized  wave  profile  corresponding 
to  (3.21),  i.e.  (^/H  as  a  function  of  ^/A,  actually  only  depends  on  two  combined 
parameters  Pi  =  PA/ cos  ^0  ^•iid  P2  =  r/P.  The  shapes  of  the  periodic  wave  given  by 
(3.21)  for  different  values  of  Pi  with  fixed  P2  are  shown  in  figure  1,  whereas  the  shapes 
of  the  periodic  wave  for  different  values  of  P2  with  the  same  Pi  are  plotted  in  figure  2. 


For  every  case,  the  corresponding  values  of  m  and  n  are  given  in  the  parentheses.  For 
fixed  P2  <  0,  both  m  and  n  increase  as  Pi  increases  (see  the  caption  of  figure  1).  On 
the  other  hand,  for  the  same  Pi,  m  increases  but  n  decreases  as  IP2I  increases  (see  the 
caption  of  figure  2).  When  r  <  0  (i.e.  P2  <  0),  as  Pi  or  IP2I  increases,  the  wave  crest 
becomes  flatter,  whereas  the  trough  becomes  steeper.  When  r  >  0,  the  situation  is  just 
opposite  because  the  crest  becomes  the  trough  and  vice  versa  under  transformation 
(C,  r)  !->■  (— (,  — r).  Two  limiting  cases:  m  — >■  1  and  m  — >  0  are  discussed  as  follows. 


Figure  2.  Profiles  of  the  periodic  wave  (3.21)  for  different  values  of  P2  with  Pi  =  5.0:  ( — 

— ),  P2  =  -0.1  (m  =  0.6799, n  =  0.3987);  ( - ),  P2  =  -0.5  {m  =  0.7229, n  =  0.2869); 

( - ),  P2  =  -2.5  (m  =  0.9523,  n  =  0.0960);  (•  •  •  •),  P2  =  -5.0  {m  =  0.9942,  n  =  0.0499). 


(a)  m  — >■  1.  In  this  limit,  the  wavelength  A  becomes  infinite  long  according  to  (3.16) 
(the  close  orbit  in  the  (x,i?)-plane  becomes  a  homoclinic  orbit,  as  shown  in  figure  Al), 
but  2K(m)/A  —>■  sec^o/\/8n.  It  can  be  readily  shown  that  sn(#,m)  — )■  tanh$, 

n(n,m)/K(m)  1/(1  —  n),  and  from  (3.19)  H  — ?•  — 4nr/(l  +  n).  In  this  situation,  the 


periodic-wave  solution  (3.21)  reduces  to  the  solitary-wave  solution: 


4rra  \  sech^$ 

1  -f-  n  /  1  —  n  tanh^  #  ’ 


$  -  KX  +  K  tan  do?  -  fiT  +  $o, 


(3.23a) 


with 

K  =  J^(2n)^/^  n  =  4K^  +  3/r  tan^  ^o-  (3.23b) 

Note  that  the  height  of  the  solitary  wave  given  by  (3.23)  increases  as  n  increases  and 
has  an  upper  bound  |2r|. 

(b)  m  — >  0.  In  this  limit,  n,H  0;  the  periodic  wave  becomes  infinitesimal.  As 


m,n  — >  0,  sn($,m)  -4  sin$,  K(m)  7r/2,  and  n(n, m)/K(m)  — >■  (1  -|-  n/2),  (3.21) 
reduces  to 


C  =  --cos 


^TT  —  — 

-^(cos  6oX  -b  sin 


Vf)  +  $0 


(3.24a) 


where 

V  =  cos  ^o(3  tan^  —  47r^  cos^  ^o/A^).  (3.24b) 


We  remark  that  (3.24)  is  the  sinusoidal- wave  solution  for  the  linearized  equation  of  (3.9).^ 
The  second  family  of  periodic  solutions  to  (3.9)  is  given  by 


C  =  r  +  m),  (3.25a) 

with 

$  =  AT  sec  6o^  +  ^o  =  K  sec  Oo  (cos  6oX  +  sin  OqY  —  Vf'j  -|-  $o,  (3.25b) 

K  =  Hf{2V2m),  (3.25c) 

V  =  cos^o  (Sr^ -b  3tan^0o  -  +  1)}  ,  (3.25d) 

where  H  is  the  wave  height,  m  (0  <  m  <  1)  is  the  modulus  of  the  snoidal  wave,  K  is 
the  A’-component  wavenumber,  V  is  the  phase  speed,  is  the  angle  of  incidence  and 
$0  is  a  phase  constant.  The  wavelength  of  the  snoidal  wave  A  is  given  by 


A  =  4K{m)cos$olK  =  8V2mK{m)cosOo/H. 


(3.26) 


Upon  using  (3.26),  (3.25a)  can  be  rewritten  as 

C  =  r  +  |/fsn(^^  +  #o).  ,  (3.27) 

Given  H,  A  and  ^o,  or  more  compactly,  given  Pi  =  HA/  cos^o,  rn  can  be  uniquely- 
determined  from  (3.26).  Note  that  the  mean  interface  corresponding  to  the  snoidal  wave 
(3.27)  is  located  at  2  =  r.  Unlike  the  wave  profiles  of  the  first  family  of  periodic  waves 
(see  figures  1  and  2),  the  profiles  of  the  second  family  of  periodic  waves  are  symmetric 
about  its  mean  interface  location.  Figure  3  shows  the  profiles  of  the  snoidal  wave  (3.27) 
for  different  values  of  Pi.  The  corresponding  values  of  m  are  given  in  the  parentheses. 
One  can  see  that  as  Pi  increases,  m  increases,  and  the  wave  crest  and  trough  both 
become  flatter  and  broader.  Two  limiting  cases  corresponding  to  m  — >  1  and  m  0  are 
discussed  below. 


Figure  3.  Profiles  of  the  periodic  wave  (3.27)  for  different  values  of  Pi  with  $0  =  — K(m): 

( - ),  Pi  =  10.0  (m  =  0.4808);  ( - ),  Pi  =  25.0  (m  =  0.8827);  (•■••),  Pi  =  50.0 

(m  =  0.9977).  For  comparison,  the  sinusoidal  wave  profile  (3.29)  is  also  plotted  in  this 
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(a)  m  — >■  1.  In  this  limit,  the  wavelength  of  the  snoidal  wave  becomes  infinite  long 
(the  close  orbit  in  the  (x,  i9)-plane  becomes  a  heteroclinic  orbit  (see  figure  A2)),  but 
4K(m)/A  — ^  H sec6o/^/8.  By  further  imposing  ^  0  as  $  — >  +oo,  we  obtain  the 
isolated  shock-like  solution  to  (3.9),  which  is 

C  =  r[l-tanh#],  ^  =  KX  +  K  UnSo?  -  nf  +  (3.28a) 


where 

K  =  |r|/\/2,  n  =  4K^  +  3K  tan"  ^o-  (3.28b) 

Since  (3.9)  is  invariant  under  transformation  C  i-)-  2r  —  C,  from  (3.28),  we  have 
C  =  r[l  -f  tanh  $],  which  is  the  isolation  shock-like  solution  satisfying  C  ^  0  as  $  ^  — oo. 
For  definiteness,  we  assume  that  C  — ?■  0  as  #  — >•  -foo. 

(b)  m  0.  In  this  limit,  H  0.  The  snoidal  wave  (3.27)  becomes  an  infinitesimal 
sinusoidal  wave: 


^  H  . 

c  =  r  -f  Y  sin 


27r 

LX 


(cos  OqX  -j-  sin  OqY  —  FT)  -1-  $o 


(3.29a) 


with 


V  =  cos  9q  (3r"  -f  3  tan"  —  47r"  cos"  ^o/A")  .  (3.29b) 

It  is  easy  to  show  that  (3.29)  indeed  is  the  sinusoidal-wave  solution  for  the  equation 
obtained  by  linearized  (3.9)  about  C  =  r. 

Note  that  the  dependent  variable  transformation  (3.8a)  does  not  include  any 
of  the  new  or  old  independent  variables.  Thus,  on  substitution  of  (3.8)  into  (3.21), 
(3.27)  and  (3.23),  (3.28),  we  obtain  the  corresponding  periodic- wave  solutions, 
solitary-wave  solutions  and  a  shock-like  solution  (satisfying  r^o  — )■  0  as  ^  ^  +oo),  to  the 
variable-coefficient  equation  (3.7).  The  phase  $  of  these  waves  can  be  expressed  as 


$  =  \/6k\^  +  Y 


(6£>i)'/"tanV<^-  Fi{q)dq  -  \ilDi/(CK) 
Jo  Jo  ^ 


+ 


Fo{q)  -  {6Dr)YH^neo  jy,is)ds  +  |  F,{s)d 


dq  ^  -h  #0.  (3.30) 
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In  the  moving  coordinates  (^jX,  Y),  at  different  X,  the  lines  of  constant  phase  remain 
straight  parallel  lines  in  the  (^,  F)-plane.  However,  their  direction  will  change  due  to 
the  contribution  from  Fi,  whose  relation  with  the  topography  is  described  by  (3.6). 
The  contribution  from  Fq  only  causes  the  lines  of  constant  phase  to  translate  and  hence 
changes  the  speed  of  the  wave  in  the  moving  coordinates. 

To  save  the  space,  only  the  explicit  expressions  for  the  family  of  the  solitary-wave 
solutions  and  the  isolated  shock-like  solution  are  given  here.  The  expression  for  the 
family  of  solitary-wave  solutions  to  (3.7)  is 


Vo  = 


F_2n 


sech^$ 


[a)uZl_3(l -|- n)J  1— ntanh^$’ 


(3.31a) 


where  #  is  given  by  (3.30)  with 

C\D.2\  (  2n 


1/2 


K 


,  f2  =  AK^  -|-  3/\  tan^  ^o- 


4/2(1  4-  n)  \  DiD-Zj 
The  expression  for  the  shock-like  solution  to  (3.7)  is 

D-2  [1  —  tanh  $] 


(3.31b) 


Vo  = 


4afj,D-3 


(3.32a)' 


where  $  again  is  given  by  (3.30)  with 


K  = 


C\D_2\ 


r,  a  =  AK^  +  ZKUrFeo. 


(3.32b) 


4fx{2DiD.zf^'' 

Note  that  when  £)_2  =  0(r  =  0),  i.e.  exactly  at  the  critical  depth  level,  we  have  only 
trivial  solution  Vo  =  0  which  vanishes  at  infinity. 

If  the  sidewalls  are  given  by 

rX 


YniX)  =  [(6Fi)'/"tan^o-Cj^Vi(s)d^ 

Yl{X)  =  F^{s)di 


<1<?  +  Ifi(0), 

d?  +  lT(0), 


(3.33a) 

(3.33b) 


any  travelling- wave  solution  C  =  C(^  +  tan^o^^  -  ooT)  will  satisfy  boundary  conditions 
(3.11).  Thus,  for  those  sidewalls  given  by  (3.33)  (note  that  they  are  parallel),  the 
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periodic-wave  solutions  we  have  found,  including  the  family  of  solitary-wave  solutions 
(3.31)  and  the  isolated  shock-like  solution  (3.32),  automatically  satisfy  boundary 
conditions  (2.12)  with  7  =  0. 

In  summary,  in  the  absence  of  rotation,  we  have  obtained  analytic  expressions  to 
describe  transformations  of  two  different  types  of  interfacial  periodic  waves  (asymmetric 
and  symmetric  about  the  mean  interface  location)  propagating  near  the  critical  depth 
level  and  over  a  weak  topography  with  linear  transverse  variation.  When  the  wavelength 
becomes  infinite,  these  expressions  reduce  to  (3.31)  and  (3.32),  which  describe  the 
propagations  of  interfacial  concave  (D_2  <  0)  or  convex  {D-2  >  0)  solitary  waves  and 
a  non-dissipative  jump  {D-2  <  0)  or  bore  {D-2  >  0)  in  the  same  setting,  respectively. 
The  presence  of  the  sidewalls  usually  will  interfere  with  the  propagations  of  these 
finite- amplitude  waves,  but  for  those  parallel  sidewalls  given  by  (3.33),  the  presence  of 
the  sidewalls  will  not  affect  the  wave  propagations.  The  most  striking  feature  is  that 
these  three-dimensional  waves  can  remain  intact  in  certain  type  of  curved  channels  with 
varying  topography;  the  only  changes  are  the  propagation  directions  and  the  phase  / 
speeds. 

Finally,  for  an  incident  Kelvin  solitary  wave  given  by 


>io(f.o,y)  = 


D-2n 

sech^  {\/6/T[e  4-  {6Dif^  tan^oK/C]}  exp(-/3y/C') 

afiD-3{l  +  n)_ 

1  -  n  tanh"*  {\/6/i:(e  +  {eDif'^  tan  OoY/C]} 

(3.34) 


propagating  in  a  varying  channel,  assuming  that  the  solution  evolved  from  this  initial 
condition  is  locally  confined  ((3.34)  satisfies  (2.18)  &t  X  =  0),  we  obtain  explicit 
expressions  for  the  first-order  and  second-order  invariants  X  (/3  =  0)  and  J  (see  (2.16) 
and  (2.17)),  which  are 


X  =  sign(D_2)-(W/C')'/" 

a 


\W-3j 


1  -  ni/2 


(/?  =  0), 


(3.35a) 


J  = 


Dr  1^-2 1 


^ZD-3  )  Q^IJ,D-3 


1  14- 

1  -1-  n  2  1  — 


Xj  [exp(-2m/C)  -exp(-2/3n/C)] . 


(3.35b) 


where  all  the  functions  are  evaluated  at  A’=0.  These  two  analytic  expressions  are  useful 
to  test  the  accuracy  of  numerical  schemes  for  (2.11)  and  (2.12). 

4.  Kelvin  solitary  wave  propagation  in  a  rotating  uniform 
channel 

Using  the  rotation-modified  KP  (rmKP)  equation,  Katsis  &  Akylas  (1987)  and 
Grimshaw  &  Tang  (1990)  investigated  the  rotation  effect  on  the  evolution  of  an  initially 
straight-crested  Kelvin  solitary  wave  in  a  uniform  channel.  They  confirmed  that  the 
rotation  gives  rise  to  a  solitary-like  wave  whose  wavefront  is  curved  back,  which  is 
in  qualitative  agreement  with  the  experiments  of  Maxworthy  (1983)  and  Renouard, 
Chabert  d’Hieres  &  Zhang  (1987).  Katsis  &  Akylas’  numerical  results  also  showed 
that  the  wave  amplitude  decays  slowly  aa  the  disturbance  propagates  downstream. 
This  indicates  that  the  solitary-like  wave  is  not  a  wave  of  permanent  form  and  the  f 
observed  attenuation  in  experiments  is  only  partly  caused  by  the  viscous  damping.  By 
stud3dng  a  coupled  set  of  evolution  equations,  which  are  asymptotically  equivalent  to  the 
rmKP  equation,  Melville  et  al.  (1989)  explained  that  the  backward  curvature  and  the 
attenuation  along  the  channel  are  caused  by  resonant  interactions  between  the  nonlinear 
Kelvin  wave  and  linear  Poincare  waves,  which  are  generated  by  resonant  forcing  of  the 
Kelvin  wave.  Similar  conclusion  was  also  reached  by  Grimshaw  &  Tang  (1990)  using  a 
small-time  asymptotic  analysis  of  the  rmKP  equation. 

In  this  section,  we  study  Kelvin  solitary  wave  propagation  neax  the  critical  depth 
level  in  a  rotating  uniform  channel  to  examine  the  effect  of  the  cubic  nonlinearity  on 
the  wave  propagation.  Following  Grimshaw  &  Tang  (1990),  we  first  carry  out  the 
small-time  asymptotic  analysis  for  (2.11)  with  constant  coefficients  and  then  numerically 
solve  the  equation.  We  shall  show  that  depending  on  the  relative  importance  of  the 
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cubic  nonlinearity  to  quadratic  nonlinearity,  an  initially  straight-crested  Kelvin  solitary 
wave  can  curve  backwards  or  forwards,  or  remain  almost  straight-crested  as  the  wave 
propagates  down  the  channel. 

To  compare  with  the  results  obtained  by  Grimshaw  &  Tang  (1990),  we  introduce 
the  following  transformation: 

10  =  X^2C(ilDifl'T,  f  y  =  CA,  (4.1) 

which  transforms  (2.11)  and  (2.12)  (with  Yr  =  =  W)  for  waves  propagating  in  a 

rotating  uniform  channel  into 

d'^U  n2  /  \ 

drde 


where 


-b  /3u  =  0,  on  A  =  0,  / 


P  =  ^^(D./3)‘'"  >  0,  1  =  W/C. 


Note  that  the  rmKP  equation  studied  by  Grimshaw  &  Tang  (1990)  did  not  include  the/ 
cubic  nonlinearity.  Equation  (4.2)  recovers  their  equation  if  P  =  0. 

It  is  easy  to  show  that  the  following  transverse  modal  functions,  satisfying  the 
boundary  conditions  (4.3),  form  a  complete  and  orthogonal  set 


</>o(A)  =  exp(-^A), 

j  \  \  X  (31  .  mr  ^  ^ 

d’n(A)  =  cos —A - sm  — A,  n  =  l,2, _ 

I  UTT  I 


(4.5a) 

(4.5b) 


The  mode  n  =  0  is  the  Kelvin  wave  mode  and  the  remaining  modes  (n  >  1)  are  Poincare 
wave  modes.  We  now  expand  the  solution  of  (4.2)  in  the  form 


u{6,  T,  A)  =  Ao{6,  r)  exp  (-/?A)  +  B{6,  r,  A), 


(4.6a) 


where 


B{e,T,\)  =  ^An{e,r)<i>r.{\) 


(4.6b) 
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represents  all  the  Poincare  modes.  Suppose  that  initially  only  the  Kelvin  wave 
component  exists,  i.e.  .B  =  0  at  r  =  0.  So  B  can  be  expected  to  remain  small  in 
sufficiently  small  ‘time’  r.  On  substitution  of  (4.6)  into  (4.2)  and  (4.3),  and  with  the 
use  of  orthogonality,  for  small  r,  the  governing  equations  for  Kelvin  wave  component 
and  Poincare  wave  component  can  be  approximated  by 
dAo 


or  a2  oO  a2  06  o9^  a2  d6\jQ  j 

f 


(4.7) 


and 


drdO 


^  80^  dX^  /?  B  +  6e 


3Be 


-0X 


AoIb - f  Be-^^^dX 

u  \  a2  Jo 

^Al  (bc-^^  -  —  I‘  Be-^^^dX 

\  02  Jo 

02 


^02 


■20X 


04  \  d'^Al 


- 


02/  80^ 


on  A  =  0, 1, 


+  0(B2),(4.8a) 


(4.8b) 


respectively,  where 


a„  =  f  e-^^^dX  =  -4  (l  - 

Jo  np  ^  ' 


(4.9) 


From  (4.8a),  it  is  clear  that  an  initially  straight-crested  Kelvin  solitary  wave  cannot 
maintain  its  original  shape,  because  it  will  generate  the  Poincare  waves  through 
resonant  forcing  (Melville  et  al.  1989),  which  will  further  interact  with  the  Kelvin  wave 
component. 

To  the  leading  order,  (4.7)  for  the  Kelvin  wave  component  becomes 


8t  02  80 

which  has  a  solitary-wave  solution: 

44o(^,r)  «  S(0,r) 


8Ao  ,  ^a3  8Al  ^a4  8Al  ,  d^Ao^^^ 


02  80 


osech^# 


1  —  T  tanh^  $  ’ 


80^ 


^  =  k{0-  4k^T), 


(4.10) 


(4.na) 
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where 


~  ~A - ^ ^  ~  ("^^3  ~  aPa^), 

Aaz  —  aPa4  2  V  2a2  ''  ^ 


(4.11b) 


and  a  (0  <  G  <  2azlPa4)  is  the  amplitude  of  the  solitary  wave.  Since  we  have  assumed 
that  5  =  0  at  r  =  0,  it  follows  from  (4.8)  that  for  small  enough  r,  the  Poincare  waves 
generated  by  the  Kelvin  solitary  wave  can  be  expressed  as 


B  =  (e-*  -  ^)]  +  O(r^).  (4.12) 


Note  that  (4.12)  satisfies  the  boundary  conditions  (4.8b)  (up  to  0(r^)),  because 


g-2/3A  _ 


€ - e 


“  =  f;«A)  ^ I' 4,HX)d\, 


(4.13a) 


(4.13b) 


Substituting  (4.12)  into  (4.7),  we  find  that 


dAo  ^azdAl  a^dAl  d^Ao  r\  2/ 

a7  +  ^7,-W  - '^7,-50  ^ -aW  - 


9P(a5  -  «3‘‘th2)-g^  +  5^^(“6  -  7ja2)-^\  =  Oij^] 


(4.14) 


The  solution  to  (4.14),  correcting  up  to  O(r^),  is  given  by 


7-2  f  Q2qZ  pa  04 

^  12(04-02/02)-^ -9^(05-03014/02)-^ 

Q  /12c>5‘ 

+  \(a,-4h2)P^-^  +0(t=)- 


(4.15) 


Substituting  (4.15)  and  (4.12)  into  (4.6a);  we  find  that  the  solution  to  (4.2)  and  (4.3), 
correcting  up  to  0(r2),  can  be  written  as 


_  osech^  {A:[^  —  Ak'^r  —  Vi{9,  A)r]} 


1  -  T  tanh2  {k[e  -  4k^T  -  Vi{e,  A)r]} 


+  o(t% 


where 


^(^.A)  =  6S(«,0)  _  3S^0,O)P  ( -  - 


(4.16) 


(4.17) 


which  can  be  interpreted  as  a  correction  to  the  incident  wave  speed  V  =  4k^  due  to 
the  Poincare  waves  (4.12).  Note  that  when  P  =  0,  the  above  results  agree  with  the 
corresponding  results  obtained  by  Grimshaw  &  Tang  (1990). 

There  are  two  terms  in  (4.17).  The  first  term  comes  from  the  quadratic  nonlinearity, 
whereas  the  second  term  comes  from  the  cubic  nonlinearity.  At  the  wall  A  =  0,  the  first 
term  is  positive,  whereas  the  second  term  is  negative  (note  that  P  >  0).  As  A  increases, 
the  first  term  decreases  and  becomes  negative  at  the  other  wall  A  =  /  (assuming  /?  >  0), 
whereas  the  second  term  behaves  just  opposite.  Therefore,  the  cubic  nonlinearity  tends 
to  cause  the  wavefront,  which  is  initially  straight-crested  in  the  transverse  direction, 
to  curve  forwards  relative  to  the  wall  where  it  initially  has  the  largest  amplitude  (here 
A  =  0  for  j3  >  0),  whereas  the  quadratic  nonlinearity  tends  to  cause  the  wavefront 
to  curve  backwards.  For  waves  propagating  near  the  critical  depth  level,  both  the 
quadratic  and  cubic  nonlinearities  are  equally  important.  The  combined  effect  on  the 
curvature  of  the  wavefront  depends  on  the  ratio  of  the  first  term  to  the  second  term  in 
(4.17).  From  (4.2)  and  (4.10),  one  can  see  that  in  the  (r,  A)  coordinates,  for  fixed  / 
incident  wave  amplitude  a,  rotation  rate  j3  and  channel  width  /,  the  relative  importance 
of  the  cubic  nonlinearity  to  the  quadratic  nonlinearity  can  be  simply  measured  by  the 
parameter  P.  Figure  4  shows  the  variation  of  the  relative  speed  correction  along  the 
crest  line  {0  =  0),  Vi(0,  A)/!^,  across  a  channel  for  a  Kelvin  solitary  wave  given  by 
(4.11)  for  different  values  of  P  with  a  =  2.0,/?  =  1.0  and  I  =  1.25.  The  corresponding 
mecin-squared  relative  speed  correction  /d[Vi(0,  A)/V]^dA  achieves  its  minimum  value 
at  Po  =  0.791  and  increases  as  |P  —  Po|  increases.  Note  that  for  the  given  a  and  (31,  there 
is  a  upper  limit  for  P,  which  is  2azl{aa4)  =  1.31.  As  P  gets  close  to  this  upper  limit,  the 
solitary  wave  (4.11)  becomes  broader  and  flatter.  From  figure  4  and  the  behavior  of  the 
mean-squared  relative  speed  correction  along  the  crest  line  as  P  varies,  one  can  expect 
that  for  small  value  of  P,  the  quadratic  nonlinearity  dominates  the  cubic  nonlinearity 
and  the  Poincare  waves  generated  by  the  Kelvin  solitary  wave  will  cause  the  wavefront 


to  curve  backwards,  whereas  for  large  value  of  P,  the  cubic  nonlinearity  dominates  the 
quadratic  nonlinearity  and  the  wavefront  will  curve  forwards.  For  some  intermediate 
value  of  P  (near  Po),  the  contributions  from  the  quadratic  and  cubic  nonlinearity' almost 
balance  each  other  and  the  Poincare  waves  generated  by  the  Kelvin  solitary  wave  are 
very  weak.  As  a  result,  the  speed  correction  due  to  the  Poincare  waves  is  not  significant, 
as  indicated  by  the  dotted  line  in  figure  4.  In  this  situation,  the  wavefront  almost 
remains  straight-crested  across  the  channel. 
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Figure  4.  A  plot  of  Vi(0,  A)/I/  as  a  function  of  A  for  different  values  of  P:  (—  •  — ), 
P  =  0.4;  (••••),  P  =  0.8;  and  (-  -  P  =  1.1. 

The  above  analysis  is  valid  only  for  small  ‘time’  r.  To  see  if  the  prediction  is  also 
applied  to  finite  r,  we  numerically  solve  (4.2)  with  (4.3).  The  numerical  scheme  developed 
by  Chen  &  Liu  (19966)  for  the  uKP  equation  is  modified  to  solve  (4.2)  with  (4.3).  The 
incident  Kelvin  solitary  wave  is  given  by  (4.11)  with  a  —  2.0,/?  =  1.0  and  I  =  1.25 
(unless  otherwise  noted),  which  satisfies  the  constraint  f^ud^  =  P(r)  exp(— /3A)  at 
T  =  0.  In  the  numerical  computations,  A0  =  0.2,  AA  =  0.05,  and  At  =  0.25  x  10“^  has 
been  used.  In  each  of  the  numerical  examples  given  below,  the  relative  errors  between 


the  numerically  evaluated  second-order  integral  invariant  and  its  exact  value  at  different 
r  are  never  in  excess  of  0.5%. 

For  small  time  r,  we  find  that  the  numerical  results  agree  very  well  with  the 
predictions  given  by  (4.16)  and  (4.17).  Figure  5  shows  the  contour  plots  of  u(r,  ^,A)  at 
finite  time  r  =  1.5  for  different  values  of  P:  0.4,  0.8  and  1.1.  The  numerical  results  show 
that  at  finite  time,  the  curvature  of  the  wavefronts  is  consistent  with  the  predictions 
based  on  the  small-time  analysis.  For  small  value  of  P  (see  figure  5a),  the  quadratic 
nonlinearity  is  dominant;  the  wave  pattern  looks  very  similar  to  that  described  by  the 
rmKP  equation  (i.e.  P  =  0).  The  curved-back  front  moves  as  a  whole  and  is  trailed 
by  a  train  of  small-amplitude  Poincare  waves,  which  is  dominated  by  its  first  mode 
(see  figure  5a).  The  maximum  wave  amplitude  attenuates  gradually  along  the  sidewall 
(A  =  0)  eis  the  disturbance  propagates  downstream  (see  dashed-dotted  line  in  figure  6). 
As  the  rotation  rate  increases,  the  extent  of  the  wavefront  curvature  and  the  amplitude 
of  the  trailing  Poincare  waves  both  increase.  For  large  value  of  P  (see  figure  5c),  the 
cubic  nonlinearity  is  dominant;  the  wavefront  curves  forwards,  accompanied  by  a  smaller. 
Poincare  wave  train  (compared  with  those  associated  with  small  P  value;  cf.  figure  5a 
with  5c).  The  maximum  wave  amplitude  along  the  channel  almost  remains  constant  (see 
dashed  line  in  figure  6).  For  P  =  0.8,  the  wavefront  is  almost  straight-crested  across  the 
channel,  with  very  small  trailing  Poincare  waves  (see  figure  56).  Note  that  the  leading 
wave  is  nearly  symmetric  about  the  crest  line.  The  maximum  wave  amplitude  along  the 
channel  decreases  at  the  beginning,  but  pretty  soon  reaches  a  constant  value  as  shown 
in  figure  6  (dotted  line). 


Figure  5.  Contour  plots  of  u(t,0,A)  at  r  =  1.5  for  different  P  values:  (a)  P 
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Figure  6.  The  variation  of  the  maximum  wave  amplitude  at  the  right-hand  sidewall 

(A  =  0)  along  the  channel  for  different  P  values:  (—  •  — ),  P  =  0.4;  ( - ),  P  =  0.8;  (-  - 

-),  F=l.l. 


Figure  7.  The  decay  of  the  amplitude  of  the  leading  wave  at  r  =  1.5  across  the  channel 

along  the  wavecrest  for  different  P  values:  (—  •  — ),  P  =  0.4;  ( - ),  P  =  0.8;  ( - ), 

P  —  1.1.  For  comparison,  the  liner  Kelvin  wave  decay  ( - )  is  also  shown. 


Figure  7  shows  the  decay  of  the  amplitude  of  the  leading  wave  at  r  =  1.5  across  the 
channel  along  the  wavecrest  for  different  P  values.  For  comparison,  the  linear  transverse 


decay  exp(— /3A)  is  also  plotted  in  the  figure  (solid  line).  In  each  case,  the  broken  line 
is  very  close  to  a  straight  line.  This  indicates  that  the  variation  of  the  wave  amplitude 
along  the  wavecrest  still  remains  exponential  for  a  Kelvin  solitary  wave  propagating 
near  the  critical  depth  level.  As  P  increases,  i.e.  the  cubic  nonlinearity  enhances,  the 
decay  rate  increases  and  approaches  to  the  linear  decay  rate.  Note  that  for  P  =  0.8 
case,  since  the  leading  wave  is  almost  straight-crested  and  its  maximum  wave  amplitude 
at  T  =  1.5  is  smaller  than  that  of  the  incident  wave  (see  figure  6),  the  decay  along  the 
wavecrest  should  be  slower  than  the  linear  decay  owing  to  the  conservation  of  energy. 
This  is  consistent  with  the  results  shown  in  figure  7  (cf.  dotted  line  with  solid  line). 


e 


Figure  8.  Contour  plot  of  u(r,  A)  at  r  =  1.5  for  P  —  1.1  in  a  wider  channel  A  =  2.5 
with  (5  =  1.0  and  a  =  2.0. 


Because  Kelvin  solitary  waves  exponentially  decay  across  the  channel,  even  if  the 
cubic  nonlinearity  is  dominant  near  the  wall  where  it  has  the  largest  wave  amplitude 
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(here  A  =  0  for  /?  >  0),  the  quadratic  nonlinearity  will  dominate  the  cubic  nonlinearity 
near  the  other  wall  if  the  channel  is  wide  enough  and/or  the  rotation  is  sufficiently  fast. 
In  this  situation,  the  wavefront  across  the  channel  will  first  curve  forwards  and  then 
curve  backwards  as  clearly  shown  in  figure  8. 

5.  Concluding  remarks 

Assuming  that  the  effects  of  nonlinearity,  dispersion,  transverse  modulation, 
rotation  and  variations  of  topography  and  sidewalls  are  weak  but  equally  important,  we 
have  derived  the  evolution  equation  for  three-dimensional  interfacial  wave  propagation 
near  the  critical  depth  level,  where  the  cubic  nonlinearity  is  comparable  to  or  dominates 
the  quadratic  nonlinearity  in  the  parametric  regime  0(e)  =  0(ju).  The  resulting 
equation  cannot  pass  the  Painleve  PDE  test  due  to  the  appearance  of  the  cubic 
nonlinear  term.  Therefore,  the  search  for  analytical  solutions  becomes  mathematically 
challenging.  Under  the  same  conditions  as  those  for  the  uKP  equation  to  be  completely 
integrable  (i.e.  no  rotation  exists,  the  variation  of  topography  is  weak  and  behaves  like  ar 
linear  function  in  the  transverse  direction),  we  have  found  two  families  of  periodic- wave 
solutions,  expressed  in  terms  of  the  snoidal  function,  to  the  variable-coefficient  evolution 
equation.  A  family  of  solitary-wave  solutions  and  an  isolated  shock-like  solution,  which 
may  be  regarded  as  a  non-dissipative  bore  or  jump,  have  also  been  obtained  «ls  the 
limiting  cases  of  the  periodic-wave  solutions.  The  analytical  solutions  show  that  these 
finite-amplitude  waves  can  remain  intact  in  certain  type  of  curved  channels  with  varying 
topography.  The  integral  invariants  associated  with  the  evolution  equation  for  interfacial 
waves  propagating  in  a  varying  channel  have  also  been  found  and  they  turn  out  to  be 
the  same  as  those  of  the  uKP  equation. 

Through  the  small-time  asymptotic  analysis  and  numerical  study,  we  have  found 
that  depending  on  the  relative  importance  of  the  cubic  nonlinearity  to  quadratic 
nonlinearity,  a  straight-crested  Kelvin  solitary  wave  propagating  in  a  rotating  uniform 


channel  will  curve  backwards  or  forwards,  trailed  by  a  train  of  small  Poincare  waves. 
When  these  two  nonlinearities  are  almost  in  balance,  the  Kelvin  wave  remains  almost 
straight-crested,  with  very  small  trailing  Poincare  waves.  The  decay  of  the  wave 
amplitude  along  the  wavecrest  across  the  channel  remains  almost  exponential  with  the 
decay  rate  increasing  as  the  cubic  nonlinearity  enhances.  The  maximum  wave  amplitude 
along  the  channel  at  the  same  location  also  increases  as  the  cubic  nonlinearity  increases. 

The  effect  of  topographic  variation  on  Kelvin  solitary  wave  propagation  near  the 
critical  depth  level  is  under  investigation  and  the  results  will  be  reported  in  the  future. 
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Appendix 

In  this  appendix,  we  shall  find  all  the  bounded  solutions  to  the  set  of  equations 

dx  „  dt? 


=  x  -px  +  q, 


dcf)  ’  d<^ 

which  is  a  conservative  system  and  has  the  following  first  integral: 


(Al) 


lixT  -  -  2px^  +  4qx)  =  lixT  -  fix)  =  C',  (A2) 

where  '  =  d/d<^,  /(x)  =  Kx^  ~  2px^  +  4qx)  and  C  is  the  constant  of  integration. 
Since  (Al)  is  invariant  under  transformation  (x^i^w)  ^  ~q)i  without  loss  of 

generality,  we  can  assume  q  >  0. 

For  the  conservative  system  (Al),  besides  the  equilibrium  solutions,  any  other 
bounded  solutions,  if  exist,  correspond  to  close  orbits  in  the  phase  plane  (X:^)-  The 
existence  of  close  orbits  depends  on  the  values  of  the  parameters  p  and  q.  Carrying 
out  analysis  in  the  phase  plane,  we  find  out  only  when  p  >  0  and  0  <  q  <  |p(p/3)^'^^, 
there  exist  close  orbits  in  the  phase  plane.  Under  these  conditions,  the  cubic  polynomial 
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•  fix)  =  “  PX  +  ?  =  0  three  distinct  real  roots,  denoted  as  q,/?  and  7  with 

7  <  0  <  /?  <  a,  and  system  (Al)  has  three  fixed  points  (a,  0),  (/?,0)  and  (7,0),  where 
(o;,  0)  and  (7,0)  are  saddle  points  and  (/?,  0)  is  a  center.  The  explicit  expressions  for  the 
0  periodic  solutions  corresponding  to  the  close  orbits  depend  on  whether  q  is  equal  to  zero 

or  not. 


Figure  Al:  The  phase  portrait  of  system  (Al)  when  p  >  0  and  0  <  q  <  |p(p/3)^'^. 


(a)  When  0  <  q  <  |p(p/3)^'^,  the  phase  portrait  is  shown  in  figure  Al.  If  the 
constant  of  integration  C  in  (A2)  satisfies  —/(a)  =  Ca  >  C  >  Cp  —  —/(/?),  there 
exists  a  close  orbit  circling  around  the  center  ((0,0).  In  this  situation,  /(x)  +  C*  =  0  has 
four  distinct  real  roots:  a,  6,  c  and  d,  which  satisfy  d<'y<c<f3<b<a<a  (see 
figure  Al).  To  obtain  the  periodic  solution  corresponding  to  close  orbit  C,  we  integrate 
the  first  integral  (A2)  along  the  close  orbit  and  obtain 

Ft, - v”  U  (A3) 

Jc  —  u){b  —  u){u  —  c)[u  —  d)\  ' 
where  (f)o  is  a  constant.  This  equation  can  be  rewritten  as: 


F{S,m)  =  ±-^{(p-  ^0), 


(A4a) 
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where 


S  =  arcsin 


{h-d){x 


[(6-c)(x-o!) 


1/2 


m  = 


(6  -  c)(a  -  d) 
[a  —  c)(6  —  d)  ’ 


(A4b) 


and  F,  the  elliptic  integral  of  the  first  kind,  is  defined  as 

du 


F(S,m)  =  [  - - 

Jo  n  — 


/I  •  2  \l/2' 

(1  —  msm  u) 

Thus,  according  to  the  definitions  of  elliptic  functions,  we  have 

(6  -  d){x  -  c) 


(A4c) 


{b-c){x-d) 


srr[K{(f)-  4>o)^m], 


where 


A^  =  [(a-c)(6-d)/8]^/". 

Hence,  the  periodic  solution  corresponding  to  close  orbit  C  can  be  written  as 

X  =  d  +  {c  — d)  -  nsn^[K{(b  —  (l>o),m]j  , 

where 


(A5a) 


(A5b) 


(A6a) 


n=^{b-c)f{b-d).  (A6b) 

Note  that  0  <  re  <  m  <  1  according  to  (A4b)  and  (A6b).  From  the  relation  between  the 
coefficients  and  the  roots  of  the  following  polynomial: 


we  have 


and 


-  2px^  +  4gx  +  4C  =  (a  -  x)(&  -  x)(x  -  c)(x  -  d), 


d  b  d  ^  c  —  0, 


(A7) 


(A8a) 


db  +  (g  -j-  h)(c  -|-  d)  T  cd  —  — 2p.  (A8b) 

Using  (A8a)  and  introducing  the  wave  weight  H  measured  vertically  from  trough  to 


crest, 


H  =  b-c, 


(A9) 
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we  can  express  a,  6,  c  and  d  in  terms  of  m,  n  and  H.  These  expressions  are  given  here: 


a  =  H  [m(l  —  2n)  +  n(2  —  n)]  [4n(m  —  n)]  \  (AlOa) 

b  =  H  [m(l  +  2n)  —  n(2  +  n)]  [4n(m  —  n)]~^ ,  (AlOb) 

c  =  H  [m(l  —  2n)  —  n(2  —  3n)]  [4n(m  —  n)]~^ ,  (AlOc) 

d  =  H  [— m(3  —  2n)  +  n(2  —  n)]  [4n(m  —  .  (AlOd) 

On  substitution  of  (AlO)  into  (A5b),  the  expression  for  K  can  be  rewritten  as 

K  =  H^{1  —  n)[8n(m  —  (All) 


The  corresponding  periodic  solution  to  (3.9),  C  =  ^  +  X  (see  (3.14b)),  can  be  obtained 
by  substituting  (AlO)  into  (A6a): 

C  =  r  —  i7[m(3  —  2n)  —  n(2  —  n)][4n(m  —  n)]“^ 

+  i7(l/n  —  1)  |l  —  nsn^[ii'(<?i>  —  ^o),”^]}  ,  </)  =  X  +  lY  —  ujT.  (A12) 

Substituting  (AlO)  and  the  expression  for  p  (see  (3.14b))  into  (A8b),  we  obtain 

a;  =  3r-2  +  -  An  +  An^) 

—  2mn(2  —  n  +  2n^)  +  n^(4  —  4n  +  3n^)][4n(m  —  n)]”^,  (A13) 

which  indeed  is  the  dispersion  relation  for  the  periodic  wave  given  by  (A12). 

(b)  When  i?  =  0,  the  phase  portrait  is  symmetric  about  the  ?9-axis  (see  figure  A2); 
/?  =  0  and  a  =  — 7  =  When  0  <  C  <  =  — /(p^^^),  there  exists  a  close  orbit 

surrounding  the  origin.  In  this  situation,  /(x)  +  C  =  0  has  four  distinct  real  roots,  i.e. 
±a,  ±6  with  b  <  a. 
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Since  the  orbit  is  symmetric  about  the  x-  ^.nd  i9-axis,  the  periodic  solution 
corresponding  to  close  orbit  C  can  be  obtained  by  integrating  the  first  integral  (A2) 
along  the  orbit  in  the  first  quadrant: 


fo  1(0=  - 


(A14) 


The  integral  on  the  left-hand  side  again  can  be  recast  in  terms  of  the  elliptic  integral  of 


the  first  kind  F: 


F(arcsm(x/6),6/a)  =  -  (M, 


(A15) 


which  can  be  rewritten  as: 


(A16) 


For  convenience,  the  constants  a  and  b  are  expressed  in  terms  of  the  wave  height  H  and 
the  modulus  m  (0  <  m  <  1)  of  the  snoidal  function: 


b  =  H/2, 


a  =  bjm  =  Hf{2m). 


(A17) 


The  corresponding  snoidal-wave  solution  to  (3.9)  becomes 


(  =  r  +  x  =  r-  +  |//sn[iJ/(2^/2m)((^-<^o)],-  <l>  =  X  +  lY  -  ujf.  (A18) 
From  the  relation  between  the  coefficients  and  the  roots  of  the  following  polynomial: 

X"  -  2px"  +  4C  =  (a^  -  xW  -  X%  (A19) 

we  have  the  following  dispersion  relation  for  the  snoidal  wave  (A  18): 

a;  =  3r2  +  3/'  -  (a^  +  6^)/2  =  +  1).  (A20) 
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